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Abstract 

We calculate the high-temperature expansion of the 2-point function up to 
order 800 in (3. We show that estimations of the critical exponent 7 based on 
asymptotic analysis are not very accurate in presence of confluent logarithmic 
singularities. Using a direct comparison between the actual series and the 
series obtained from a parametrization of the form (/3 C — /3)~ 7 (Ln(/? c — /?)) p +r), 
we show that the errors are minimized for 7 = 0.9997 and p = 0.3351, in very 
good agreement with field-theoretical calculations. We briefly discuss the 
related questions of triviality and hyperscaling. 
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I. INTRODUCTION 



The dimension four plays a doubly important role in physics. First, it is the dimension 
of space-time which is relevant for a relativistic description of a large class of phenomena, 
from electricity and magnetism to scattering processes at the highest experimentally ac- 
cessible energies. Second, it is the upper critical dimension for scalar field theory. If one 
analytically continues the renormalization group equations (usually derived within some 
approximation) to non-integer dimensions, it appears that when the dimension tends to 
four from below, the non-trivial fixed point merges with the Gaussian one. This justifies the 
e-expansion. 

It is thus commonly accepted that in four dimensions, the critical exponents are the 
trivial ones (i.e. those obtained from mean field). Unfortunately, it often difficult to find 
clear evidence for or against trivial exponents, for instance, from high-temperature (HT) 
series [||,[3J or a finite volume calculation [|J]. The root of the problem is the existence of 
a marginal direction which makes the approach to the fixed point more intricated than in 
three dimensions. The corrections to the power laws can in principle be obtained from the 
Callan-Symanzik equations, provided we know the exact form of the various functions (beta, 
gamma, ... ) entering into them. Using the lowest order in perturbation theory, Brezin, 
Le Guillou, and Zinn- Justin || found that the trivial power divergences get multiplied by 
rational powers of Ln((3 c — (3). It is important to check this result with methods independent 
of perturbation theory. In particular, it is conceivable that there exist non-trivial fixed points 
which cannot be revealed by perturbation theory. 

The technical challenge which appears in any kind of calculation is to distinguish between 
a small change (with respect to the trivial value) in the critical exponent and a slowly varying 
(compared to the trivial singularity) multiplicative change. This difficulty appears clearly in 
the asymptotic analysis of the high-temperature expansion of the susceptibility, where the 
leading term of the extrapolated slope defined in Eq. Q3.4|) (7 — 1) can be small compared 
to corrections proportional to the inverse of the logarithm of the order, unless one can reach 
an astronomically large order. 

Another interesting feature of the field-theoretical method is the so-called hyperscaling 
relation among the power singularities of the 2- and 4-point (subtracted) Green's functions 
at zero momentum. In three dimensions, the violations of hyperscaling |J are hard to resolve 
by high-temperature calculation. This is still a controversial topic. In four dimensions, 
conflicting |||3| conclusions were drawn from the high-temperature series. 

The confirmation of the field-theoretical results would require that an unbiased estimate 
of the main power singularity and the power of the logarithmic correction come close to their 
predicted values, with errors compatible with (small) higher-order corrections. We propose 
here to test the field-theoretical results using an expansion in the kinetic term (also called 
high-temperature expansion), in a model which is obviously non-trivial in three dimensions, 
but where calculations are easier than in nearest-neighbor lattice models. 

The hierarchical model j8| is a non-trivial approximation of models with short range in- 
teractions, which is well-studied [§,0 , and for which we can calculate the high-temperature 
expansion Jill to a very large order. The recursion relation which summarizes the renormal- 
ization group transformation of this model is closely related to the approximate recursion 
formula discussed by Wilson The qualitative and quantitative aspects of this relationship 
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are discussed in Ref. |L2] . 

In recent publications [13,14 , we reported results concerning the high-temperature 
expansion of Dyson's hierarchical model in three dimensions. We calculated the high- 
temperature expansion of the magnetic susceptibility up to order 800 with Ising and Landau- 
Ginzburg measures. This allowed us to obtain a value of the critical exponent 7 of 1.300 
in D = 3, with estimated errors of order 0.002. This result is consistent with the results 
obtained with the e-expansion @.[K]. 

We found clear evidence for oscillations in the quantity, called the extrapolated slope 
|T5[ (see section below), used to estimate the critical exponent 7. When using a log scale 
for the order in the high-temperature expansion, these oscillations become regularly spaced. 
Our interpretation of the data was consistent with the hypothesis that the eigenvalues of the 
linearized renormalization group transformation are real, but that the constants appearing 
in the conventional parametrization of the magnetic susceptibility should be replaced by 
functions of j3 c — (3 invariant under the rescaling of [3 C — (3 by Ai, the largest eigenvalue of 
the linearized renormalization group transformation. This possibility has been mentioned in 
the past by K. Wilson [|I|] and developed systematically by Niemeijer and van Leeuwen |L6| . 
Our analysis provided good evidence that the oscillations appear with a universal frequency 
in good agreement with theoretical expectations, but with a measure-dependent phase and 
amplitude. 

Subsequently, more efficient methods of calculation, based on finite dimensional projec- 
tions of the Fourier transform of the recursion formula, were developed. As explained in 
detail in Ref. |17]], the effects of such truncations can be controlled with a precision which 
is better than exponential when the dimension of the truncated space increases. 

In this paper, we study the high-temperature expansion of Dyson's hierarchical model 
in dimension 4. For the sake of completeness, we briefly review the method of calculation in 
section ||. The conventional methods |15| , |18| used to estimate the critical temperature and 
a critical exponent from a high-temperature series are reviewed in section fTJ. We show that 
in the presence of logarithmic corrections to the scaling laws, the asymptotic behavior of the 
corrections is modified. The extrapolated ratio defined in Eq. (|3.3p provides an estimate of 
the critical temperature with corrections of order rrT 1 x (Ln(m))~ 2 , where m is the order in 
the high-temperature expansion. In the following, we continue to use the notation m with 
the same meaning. Using the expansion of the susceptibility up to order 800, we obtained a 
value of the critical temperature which agreed with the high-precision determination of Ref. 



T7J with errors of less than one part in 10,000. On the other hand, the extrapolated slope 



defined in Eq. (g) estimates the critical exponent minus one with corrections which are only 
suppressed by (Ln(m)) -1 . If this weak suppression is not recognized, one may conclude that 
the critical exponent 7 takes a value larger than the trivial one. More generally, asymptotic 
analysis is not adequate to distinguish between a value of 7 close to 1 and a correction to 
the scaling laws which is less singular than a power. 

In section |iy], we analyze the high-temperature expansion of the susceptibility without 
relying on the asymptotic behavior of the coefficients. We use h(m) = {r m /3 c — l)m, a 
function which represents the difference between the ratio of successive coefficients r m and 
its asymptotic value . The function h(m) can be calculated exactly using either the 
empirical series or the series corresponding to a given assumption on the analytical form of 
the susceptibility. Taking the sum over a large range of m of the square of the differences 
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between these two values of (h(m)) 1 , one can get an error function which indicates how 
good the analytical assumption is. We found that the parametrization 

X = (& - /3)^(A (-/3- 1 Ln(l - + At) (1.1) 

provides very good fits of the data for 7 ~ 1 and p ~ |, which is the field-theoretical f| 
result. In order to decide how accurate the agreement is, we have considered fixed values 
of 7 in the vicinity of 1 and equally spaced by 1CT 4 steps. For each of these values, we 
have determined the values of p and Ai/Aq which minimize the error function. This error 
function behaves like a paraboloid near its minimum at 7 = 0.9997 and p = 0.3351, in good 
agreement with the field-theoretical calculation. The errors on this estimate are mostly 
systematic. To get more accurate results, one needs to replace the constant A\ by a slowly 
varying function. 

Another quantity which can be studied using the high-temperature expansion is the 
dimensionless renormalized coupling constant [|19|], denoted A4 hereafter, obtained by multi- 



plying the connected four-point function at zero momentum by the eighth (D + 4 in general) 
power of the renormalized mass. For D < 4, this quantity is designed to have a finite and 



non-zero limit when (3 — > (3 C . In the case D = 4, we have checked with good accuracy (17 
that A4 goes to zero like (Ln(f3 c — /?)) _1 for the model studied here. The calculation of the 
HT coefficients of A4 involves the subtraction of the disconnected part and suffers the same 



type of numerical problems as the direct calculation of A4, as discussed in Ref. fl7|| . For this 
reason, we were only able to extract a series of 30 coefficients. The analysis of this series is 
consistent with the fact that A4 goes to zero when (3 — > j3 c (triviality), but it is not possible 
to distinguish a (Ln((3 c — approach to zero from a {(3 C — (3) 1 ^ 2 approach, which would 

be necessary to establish whether or not hyperscaling holds. This question has been settled 



in Ref. [ 171 , and this section illustrates the inconclusiveness of results obtained from short 
series. 

In conclusion, we have shown that by using sufficiently long series and methods of analysis 
not relying on an asymptotic expansion, it is possible to obtain very good agreement between 
calculations based on field theory and those based on high-temperature expansion in the 
upper critical dimension. We emphasize that the main interest of the high-temperature 
expansion is to allow us to probe global features of the renormalization group flows which 
cannot be approached using renormalized perturbation theory or an analysis of the linearized 
behavior near the fixed point. An example of such a global feature is the existence of log- 



periodic oscillations |13,14|) which play an important role in D = 3, but have an almost 



negligible effect in D = 4, as shown in section |TJ. Another example of a global feature 



could be the existence of a non-trivial fixed point. The good agreement found in section [TV 
makes this possibility very implausible for the model studied here. 

II. CALCULATIONS OF THE HT COEFFICIENTS 

The calculation of the high-temperature expansion of the unsubtracted 2fc-point functions 
of Dyson's hierarchical model can be performed iteratively using the basic recursion formula 



in its Fourier form |TTJ. This method has been discussed extensively in Refs. [|13| , |T4|| . For 
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the sake of being self-contained, we briefly explain the basic method of calculation. More 
details, justifications, and motivations can be found in Refs. |pj [l3| , |T^| . 

The recursion formula for the rescaled Fourier transform R n (k) of the local measure for 
blocks of 2 n sites reads 

R n+ i(k) = C n+1 eM-\^\s 2 ) n+1 ^){R n ^)f , (2.1) 

where c is an adjustable parameter which takes the value 2 1_ "E>, in order to approximate 
.D-dimensional models. In the following, we will only consider the case D = 4, which means 
c = a/2. The rescaling operation commutes with iterative integrations, and the rescaling 
factor s can be fixed at our convenience. In order to obtain stabilized expressions in the 
high-temperature phase, we will take s = v2 in the following. We fix the normalization 
constant C n in such way that Rn{0) = 1- R n {k) then has a direct probabilistic interpretation. 
If we call M n the total field 4>x inside blocks of side 2 n , and < ... > n the average calculated 
without taking into account the interactions among these blocks, we can write 

Uk) = ±t^L^^. ( , 2) 

We see that the Fourier transform of the local measure obtained after n iterations generates 
the zero-momentum Green's functions calculated with 2™ sites. All the calculations done 
here use an initial Ising measure, which means that Ro{k) = cos(k). Since we are interested 
in the leading singularity, this choice should play no role |19| in the discussion. 

The high-temperature expansion of the zero-momentum Green's function can be obtained 
from an expansion of Eq. ( |2.1|) in powers of (3. The most important sources of errors are 
the round-off errors. After 100 iterations, the relative errors on the mth coefficient |L4] are 
of the order of m x 10 -15 . With the choice s = \/2, the coefficients reach a finite value in 
the infinite volume limit. Actual computations are made at large but finite volume (i. e. at 
finite n). The relative difference between the coefficients at finite and infinite n goes to zero 



ll| like (|) n . For D = 4, the choice n = 100 means that (|) n = 2 50 , which is smaller than 



the numerical errors. 

Such a calculation is in general time-consuming when one wants to calculate more than 
100 coefficients. It is, however, possible to save time by using finite dimensional approxima- 



tions [IT] of degree I for the generating function: 

R n (k) = 1 + a nA k 2 + a n , 2 k 4 + + a n jk 21 , (2.3) 



with / much smaller than the required dimension m+1 necessary for an exact []TTj] calculation. 
After each iteration, non-zero coefficients of higher order (a n+ i^ + i etc. ) are obtained, but 
set to zero in the next iteration. The /-dependence of the high-temperature coefficients of 
the susceptibility is discussed in Ref . JL7) . If b$ denotes the value of b m in a truncated space 
of dimension I, we found that 

where s and % are, respectively, the slope and intercept of the corresponding fitted line, 
as shown in Fig. |l|. The intercepts are approximately 2.3, while the slopes depend on m. 
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Eq. (j2.4|) represents suppressions which are better than exponential. From this figure, we 



can check, for instance, that for m = 400 (which is the maximal value used in section IV). 
the extrapolated errors at Z = 40 are significantly lower than the numerical errors. Using 
extrapolation in m, it was estimated in Ref. |L7j that in the case D = 4, I = 38 was sufficient 
to calculate 61000 • 

In summary, the following calculations will be performed with I = 50 and n = 100. The 
above discussion shows that this choice guarantees that the systematic errors are smaller 
than the numerical errors. 



III. THE LIMITATION OF THE ASYMPTOTIC ANALYSIS IN PRESENCE OF 
CONFLUENT LOGARITHMIC SINGULARITIES 

In this section, we study the singularities of the susceptibility using its high-temperature 
expansion 

00 

x(P) = E (3.1) 

m=0 

We define r m = 6 m /fe m _i, the ratio of two successive coefficients. When D < 4, one expects 
Ij| that 



X = (ft - P)- J (A + A 1 {f3 c - P) A + ....) , (3.2) 



and it is convenient to introduce quantities [15] called the extrapolated ratio (R m ) and the 



extrapolated slope (S m ) in order to estimate (3 C and 7. These quantities are defined as 

R m = mr m - (m - l)r m _i , (3.3) 

and 

S m = mS m - (m - l)S m -! , (3.4) 

where 

S m = -m(m - l)(r m - r m _i)/(mr m - (m - l)r m _i) (3.5) 

is called the normalized slope. When A$ and A\ are constant, one finds [15| that the 1/m 
corrections disappear: 

S m = 7 -1- Bm- A + 0(m- 2 ). (3.6) 



However, for the hierarchical model in D = 3, large oscillations were observed [13 
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S m and it was recognized [|T3| , p!4| that A and Ai should be considered as functions of 
(3 C — (3 invariant under the rescaling of (3 C — (3 by Ai, the largest eigenvalue of the linearized 
renormalization group transformation. The asymptotic analysis (when m becomes large) of 



the extrapolated slope in this modified situation is given in section 3 of Ref. It was 



found that 1/m corrections with rather large coefficients reappeared. Nevertheless, it was 
possible to extract the critical exponent 7 with estimated errors of 0.2 percent. 
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The situation is very different in D = 4, as shown in Fig. |2|. The oscillations are barely 
visible for low values of m, and not visible at all for larger m, where S m appears to decay 
smoothly. If the parametrization of Eq. ( |3.2| ) and its corollary Eq. ( |3.6| ) applied, one might 
conclude that 7 is close to 1.05. However, if we plot the inverse of (Sm)' 1 versus Ln(m), 
we find the linear behavior shown in Fig. [| This shows that S m decays like 1/Ln(m), so 
Eq. Q3.6|) does not provide an adequate description of the situation. The deviation from the 
linear behavior shows an interesting fine structure shown in Fig. [|. For m near 400 (Ln(m) 
near 6), ones sees that the amplitude of oscillation is almost four orders of magnitude smaller 
than (Sm)' 1 itself. For such a values of m, the numerical errors become comparable with 
the oscillations. For larger values of m, the numerical errors become larger and wash out the 
oscillations. The numerical errors on (Sm)" 1 in D = 4 are of the same order of magnitude as 
what we would estimate in D = 3 from the error analysis of Ref. |TJ] . The main difference 



is that the oscillations have a much smaller amplitude in D = 4. In the following, we will 
treat the oscillations on the same footing as the numerical errors, which is justified for m 
sufficiently large. 

We will now revisit the asymptotic analysis of R m and S m in a more general case than 
Eq. ( |3.2|) with Aq and A\ constant. Our main assumption will be that 

X((3) = (1 - jyG(l - |) , (3.7) 

where G is such that 

/Zmm -°°^P) = ° • (3 ' 8) 

This restriction includes the case where G(l — -§-) grows like a positive power of a logarithm 
when p goes to (3 C . We then proceed as in ref. |1| and explain the principle of the asymptotic 
expansion. We use the residue theorem in the complex (3 plane to get an integral representa- 
tion of the coefficients. Next we treat the integral with the steepest descent method. Using 
an exponential parametrization for the integrand of the mth coefficient, one finds that the 
phase has a maximum for a value of y = 4- such that 

Pc 

7 G'(l -y) 

^~W^v)^ m + 1 - (3 ' 9) 

The basic principle of calculation is that the second term of the 1. h. s. of this equation can 
be treated as a perturbation, for m sufficiently large. Neglecting this second term, we get 

y = i-- + o(\) . (3.10) 

Eq. ( p.8[ ) is then seen as the condition which allows us to treat the second term of the 1. 
h. s. of Eq. ( |3.10|) as a perturbation. Finally, one finds [0 that for large m, the leading 



contribution to the mth coefficient has the form 



b m oc /5- m m 7 - 1 G(— ) . (3.11) 
m 
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Before going further, we introduce a parametrization of the ratio of successive coefficients: 

r m = (h(l + -h(m)) . (3.12) 

Pc m 



This definition is independent of any kind of expansion. From Eq. |3.11| , we found the 
asymptotic estimate 

h(m) = ( 7 -l---^) + ... (3.13) 



If we consider the case 



we obtain 



and 



G{x) oc (-Ln(x)) p , (3.14) 



^ = (> + °(^W )) ' (3 ' 15) 



S m = 7 - 1 + 0(-— -) . (3.16) 
Ln[m) 



From this, we can conclude that under the assumption of Eq. ( [3.14D , asymptotic analysis 



justifies using R m as an estimator for with estimated errors on the order of 10 4 . This 

_ Pc 

quantity is displayed in Fig. |5|. As expected, no oscillations are visible. The change between 
m = 200 and m = 800 is less than 10~ 4 , which is consistent with tT \ m corrections. If 

we use -Rsoo as our best estimate, we obtain (3 C = 0.665548, which is in good agreement with 



our accurate calculation fll7[ , where we found f3 c = 0.6654955715318593. The discrepancy 
has the same order of magnitude as the small variations noted above. 

On the other hand, for S m , the corrections to 7 — 1 are not very small. For instance, 
for m = 800, (Ln(m)) -1 ~ 0.15, and it seems implausible that one could establish that 
I7 — 1| < 10~ 3 on the grounds of an expansion in this not- very small parameter. More 
generally, it takes exponentially large m for the "corrections" in (Ln(m)) -1 to become 
smaller than the "leading" 7 — 1 when this quantity is small. Thus it seems desirable 
to use non-asymptotic methods, the subject of the next section. 



IV. A DIRECT ESTIMATION OF THE CRITICAL EXPONENTS 

In this section, we propose to use direct calculations of h(m), defined as 

h(m) = (r m (3 c - l)m . (4.1) 

This quantity can be calculated exactly under some assumption regarding the susceptibility, 
and calculated exactly from the empirical series. We emphasize that h(m) is defined from 



Eq. ( 3.12 ) and its calculation does not require any kind of expansion. However, we need to 
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provide an estimate of j3 c . In the following, we will take the most accurate value [17) of (5 C 
quoted in the previous section rather than the approximate values obtained from R m . 

A simple assumption on the leading singularity of the susceptibility in D = 4 is given by 
the result of a field-theoretical calculation || 

X=(Pc-(3)~ 1 (A (\ln(P c -(3)\)^ . (4.2) 

This lowest-order result would also be obtained for Dyson's model, because at this order, 
numerical factors (integrations over the angles), which are model-dependent, cancel. Given 
that r m is the ratio of two successive coefficients, it is independent of Aq and it transform 
as r m — > r m s~ x under a rescaling (3 — > s(3. Consequently, r m (3 c is independent of the choice 
of A and [3 C . We have thus calculated h{m) from the expansion of 

f{x) = {l-x)-\- 1 ^—^) 1 * (4.3) 

x 

in x, about x = 0. The variable x stands for j3/(3 c . The division by x does not change 
the leading singularity || when x — > 1 while providing a regular expansion around x = 0. 
Under the assumption of Eq.(O), we find from Eq.( |3.1ij ) that asymptotically h(m) tends 



to a small and possibly zero constant plus a correction which decays like Ln ^ ■ It is thus 
natural to plot (/i(m)) -1 versus Ln{m). Such a plot is provided in Fig. || where we compare 
with (/z(m)) -1 calculated directly from the D = 4 HT series, using the definition Eq. ( |4.1| ). 
The two (approximate) lines are separated by an almost constant gap. We tried to modify 
the assumption Eq. ([4.3|) in such way that the two lines coincide. The only satisfactory 
solution we found was the modified assumption 

fix) = (1 - x)-\(- ln{1 ~ X) )t + r) , (4.4) 

x 

where r has to be determined by an error-minimization procedure which we now proceed to 
explain. 

For notational purposes, we call t(m) the "true" value of (/i(m)) -1 obtained from the HT 
series and a(m) the value of (/i(m)) -1 corresponding to an assumed series such as the one 
obtained from Eq. ( |4.4| ). In practical calculations, we have used the instruction Series in 
Mathematica to calculate a(m). It should be noted, that for large orders, rational values of 
the exponents give better numerical results. In addition, if the denominator of this rational 
exponent gets too large (typically 10 7 for a calculation up to order 400), one runs out of 
memory. This procedure is quite time-consuming when one goes to large order. Since such 
a calculation will have to be repeated many times in the rest of this section, we have used 
the region 300 < m < 400 to evaluate the discrepancy between a(m) and t(m). As we can 
see from the discussion of section |T|, in this range the oscillations are already small and 
the numerical errors not too large yet (see Fig. |j). We have thus determined the parameter 
r in Eq. ( |4.4j ) by minimizing 

m=400 

e = W m ) - a ( m )) 2 • ( 4 - 5 ) 

m=300 

The values of E for values of r separated by 0.001 are shown in Fig. |7|. The curve can be 
fitted very well by a parabola. The minimum of this parabola is then determined analytically 
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from the three values defining the fitting parabola. This allows us to find the value of r 
with a precision of 10~ 6 , three orders of magnitude smaller than the original resolution. The 
value of r minimizing E found with this procedure is -0.435622, corresponding to a value 
of E of order 10~ 4 . Subsequently, we checked this answer by repeating the calculation of E 
with a resolution 10~ 6 in r and found the same answer. For this value of r, t(m) and a(m) 
cannot be distinguished in a graph like Fig. |6|. The difference between a(m) and t(m) is 
shown in Fig. In the region where the fit was performed (300 < m < 400), the differences 
are 4 orders of magnitude smaller than the values themselves. 

We conclude from this analysis that Eq. Q4.4| ) is a very good guess concerning the leading 
and subleading singularities of the susceptibility. However, we would like to see if it remains 
the best guess if we allow the exponents to change. In other words, we would to see if 
different values of the exponents could also be acceptable from the point of view of the 
high-temperature expansion. We have thus considered a more general assumption: 

f(x) = (1 - x)-^((- ln{1 ~ x) ) p + r) (4.6) 

x 

and studied E function of 7, p, and r. 

Near a minimum, E behaves generically as a three-dimensional paraboloid. In this 
region, one can "eliminate" r by fixing its value in such a way that E is minimized with 
7 and p kept constant. The variable r is thus replaced by a linear combination of 7 and 
p plus a constant and we can then work with a two-dimensional paraboloid. A section of 
this paraboloid defined by the condition 7 = 1 is shown in Fig. |9|. In a second step, one 
can similarly eliminate p with 7 fixed by requiring that it takes the 7-dependent value that 
minimize E. In the case 7 = 1 illustrated in Fig. ^|, this value of p is 0.32775, not far from the 
expected |J value |. This show incidentally that a biased estimate of p is in good agreement 
with the field-theoretical result. Taking values of 7 separated by 10 -4 , we have similarly 
calculated the value of p given by the minimization condition. The results are shown in Fig. 



10] . The linear behavior was expected: since near the minimum E is a quadratic form, the 
minimization condition is linear. Using this linear relation to eliminate p, E{pf) becomes a 
parabola. The minimum value taken by this function is then the minimum of the initial 



function E("y,p,r). This function is shown in Fig. |TT]. E is minimized for 7 = 0.9997, which 
according to Fig. |TD| corresponds to a value of p of 0.3351. 

In practical calculations, it is convenient to replace parabolic fits by successive appli- 
cations of Newton's method. This method has an adjustable resolution and it allows one 
to start in regions away from the minimum, and where the parabolic behavior does not 
necessarily hold. 

It is difficult to estimate the errors on our result. Since the parabolas shown above are 
reasonably smooth, it seems unlikely that the numerical errors or the oscillations, which 
should have about the same size, play any significant role. Most likely, the main source 
of error is that r has been considered as a constant. If instead we allow r to be a slowly 
varying function of /3, we expect in a model independent way, that these slow variations in 
/3 will induce slow variations in m of the quantity (h(m))~ l . In the interval of m considered 
for the calculation of E, the slow variations can be approximated with polynomials. In 
order to get an idea about how low E could become under such a circumstances, we have 
fitted the differences between t{m) and a{m) displayed in Fig. R] and calculated the value 



10 



of E obtained after subtracting these fits from the original differences. For a linear fit, we 
obtained E = 6 x 10~ 8 and for a quadratic fit E = 8 x 10 -10 . This shows that by keeping 
7 = 1 and p — 1/3 and allowing r to to be a slowly changing function written in terms 
of parameters which are adjusted to minimize E, we can obtain values of E comparable to 
those obtained by keeping r as a constant and allowing 7, p and r to be adjusted in order 
to minimize E. For definiteness, with 7 = 1 and p and r varied to minimize E, we obtain 
E = 6 x 10 -9 . Varying 7, p and r, we obtain E = 4 x 10~ 10 . In conclusion, if we want a 
more precise estimation of the critical exponents, we also need more information regarding 
the /9-dependence of the subleading singularities. 



V. TRIVIALITY AND HYPERSCALING 

Another quantity which can be studied using the high-temperature expansion is the 
dimensionless renormalized coupling constant |T9 



A 4 = -G\m D R ^ , (5.1) 

where G\ is the the zero-momentum connected Green's function and m# the renormalized 
mass. For D < 4, this quantity is designed to have a finite limit when (3 — > /3 C . In the case 
D = 3, we have checked ]17[ by a direct calcualtion that A 4 reaches the value 1.92786 when 
(3 — > f3 c . In the case D = 4, we have checked with a good accuracy that, in the same limit, 
A4 goes to zero like (Ln(j3 c — f3))~ l . Thus we have direct evidence that in these two cases 
the power singularities cancel in Eq. ( |5.1| ) — in other words, that hyperscaling holds. 

Bearing in mind that there is no wave function renormalization (rj = 0) in the hierarchical 
model, we will define as in Ref. W7\ that A 4 is the limit where n — > 00 of 



-.: Mj > n -3(< Ml > n f 
2"( ; " )^ +2 



with the same notation as in Eq. (|2.2j ). Equivalently, with the rescaling factor fixed to 



s = \[2 in Eq. ( |2.1|) and the convention of Eq. (|2.3|) , 

A 4n = 12 Q "' 1 ~ 2a D "' 2 2". (5.3) 
(~2a nil )T+ 2 1 J 

The calculation of the HT coefficients of A4 involves the subtraction of the disconnected 
part and it suffers the same type of numerical problems as the direct calculation of A4, as 
discussed in Ref. |17j . For this reason, we were only able to extract a series of 30 coefficients. 



The quantity him) defined in Eq. (|4.1| ) corresponding to this series is displayed in Fig. |l"2| . 
The figure indicates that this quantity has damped oscillations. The average value of h(m) 
in the displayed interval is -1.4. From Eq. (|3.13|) , this is consistent with the fact that A4 has 
a finite limit when j3 — > (3 C , but it not possible to distinguish a {Ln{(3 c — approach to 

zero from a (/3 C — P) 1 ^ 2 approach. For comparison, we have displayed in Fig. |12| the function 
h(m) corresponding to the series generated by —x/Ln(l — x) and (1 — x being short for 
-§-. It is clear that the oscillations make the discrimination between these two asumptions 
impossible. 
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Another way of seeing that the series is too short to describe the details of the behavior 



near j3 c is to plot the truncated expansion of A4 up to order 30. This is done in Fig. [13 . 



The HT expansion indicates correctly that A4 goes to zero when j3 increases. However, the 



behavior near j3 c is not accurate. For comparison, Fig. [13] also shows the leading critical 
behavior estimated in Ref. [[T7|], namely 

(5.4) 



-1.955 - 0.746 x Ln(p c - p) ' 

The data interpolates nicely between the two types of behavior, but we see that there is no 
region in the figure where they overlap. The order 30 HT expansion gives accurate results 
for p c — p > 3 x 10 -2 , while Eq. fl5.4|) becomes accurate when p c — P < 10~ 3 . 

In summary, the truncated expansion makes clear that A4 goes to zero when P —>■ ft c . In 
other words, the theory is trivial. However, the series is too short to extract accurately the 
precise way it approaches zero, and one cannot decide from this information whether or not 
hyperscaling holds. 



VI. CONCLUSIONS 

There have been questions || in the past regarding possible discrepancies between field- 
theoretical calculations based on the renormalization group approach and calculations based 
on the high-temperature expansion. Using a scalar model in the upper critical dimension, 
where all the conventional expansions can be compared with direct calculations, we claim 
that the field-theoretical result concerning the leading singularity of the two-point function 
at zero momentum given in Eq. ( |4.2| ) can be reproduced very well by the high-temperature 
expansion. 

Using a parametrization of the subleading singularities depending on a single constant 
r, we obtained an optimal agreement for the choice 7 = 0.9997 and p = 0.3351. With this 



choice, the error on (/i(m)) -1 defined in section [V, is less than one part in a million for 
300 < m < 400. The small discrepancies between our estimate of the critical exponents and 
the field-theoretical values 7 = 1 and p = 1/3 are not significant because it is possible to 
show that small changes in the exponents and allowing r to slowly vary have comparable 
effects for the quality of the fit. 

The present study shows that the use of asymptotic analysis or the use of a short series 
can be misleading. Given the length of the series available, asymptotic analysis may be useful 
for order of magnitude estimates but not for an accurate determination of the exponents. 

There is still room for improvement. One could use calculations at fixed P to study the 
corrections to the parametrization of Eq. Q4.6|) . This procedure could be pursued up to the 
point where the main source of errors would be the numerical errors on the coefficient. 

The use of the high-temperature expansion allows us to probe global features of the 
renormalization group flows which cannot be approached using renormalized perturbation 
theory or an analysis of the linearized behavior near the fixed point. In particular, our 
analysis makes implausible, for the model considered here, unconventional possibilities such 
as the existence in the upper critical dimension of a non-trivial fixed point characterized by 
non-trivial exponents. 
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FIGURES 

1 — dependence of b (D = 4) 
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FIG. 1. /-dependence of the high-temperature coefficients bm calculated in truncated spaces of 
dimension I. 
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Inverse Extrapolated Slope (Ising = 4) 
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FIG. 4. Difference between 5 m 1 and a linear fit in Ln(m). 
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FIG. 5. The extrapolated ratio R m versus the order m in the HT expansion. 
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FIG. 7. The function E defined in Eq. ( [4. 5|) versus r, with a(m) calculated from Eq 
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FIG. 8. a(m) — t{m) versus m for a(m) calculated with r = —0.435622 in Eq.( |4.4| 
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FIG. 9. The function E versus p with 7=1 and r chosen to minimize E. 
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FIG. 10. Values of p minimizing E at given 7 and optimum r. 
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FIG. 11. The function E versus 7 with values of p and r chosen to minimize E. 
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FIG. 12. The function h(m) corresponding to the HT expansion of lambda^ (dots) compared 
to the same function for —x/Ln(l — x) (circles) and (1 — x)a (squares). 
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FIG. 13. A4 versus f3, exact (dots), with the HT expansion up to order 30 (dashed) and in 
leading singularity (line). 
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